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PREFACE 


This  program  was  conducted  for  the  Defense  Nuclear  Agency 
under  contract  DNA001-77-C-0280  to  provide  experimental  data  on 
shock  attenuation  in  the  M-X  trench.  The  effort  was  sponsored  by 
the  Strategic  Structures  Division  of  the  Shock  Physics  Directorate 
and  was  performed  during  the  period  5  July  1977  to  31  December 
1978.  Dr.  George  Ullrich  was  the  DNA  Technical  Representative. 

The  SAI  Program  Manager  was  Duane  Hove.  Dr.  James  Craig 
conducted  the  shock  tube  tests  and  performed  the  data  analysis. 

Dr.  Raad  Issa  developed  the  one-dimensional  computer  code  used 
for  test  predictions,  and  Dr.  Kenneth  Haines  directed  the  optical 
measurements  team. 

We  are  grateful  to  Dr.  Robert  Dannenberg  of  the  NASA 
Ames  Research  Center  for  continuing  cooperation  and  suggestions. 
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NOMENCLATURE 


A  Sound  speed 

C_  Shock  velocity 

b 

Friction  coefficient 
D  Tube  diameter 
E  Energy 

k  Roughness  height 

L  Length 

m  Molecular  weight 

M  Mach  number 

P  Pressure 

Pr  Prandtl  number 

R  Tube  radius 

R  Dimensionless  velocity  at  roughness  height 
ReD  Reynolds  number  based  on  diameter 
S  Spacing 

St  Stanton  number 

T  Temperature 

U  Velocity 

V  Volume 

X  Distance 

Y  Ratio  of  specific  heats 

r  Acoustic  impedance 

v  Kinematic  viscosity 

p  Fluid  density 

x  Shear  stress 
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1.0 


INTRODUCTION 


As  an  alternative  to  fixed  based  concepts,  the  USAF  is  con¬ 
sidering  locating  an  advanced  missile  (M-X)  in  multiple  trench  bases 
to  increase  the  ability  of  the  force  to  survive  preemptive  attack. 
Success  of  this  approach  relies  on  designing  a  cost  effective  bas¬ 
ing  system  with  sufficient  hardness  to  nuclear  attack  such  that  an 
agressor  would  have  to  exhaust  his  arsenal  to  achieve  assured 
destruction. 

Airblast  from  a  nuclear  weapon  presents  a  critical  design 
environment  for  the  M-X  trench.  Blast  waves  will  propagate  in  the 
trench  due  to  loading  during  the  cratering  process,  due  to  loading 
of  a  previously  breached  trench  or  due  to  a  burst  in  the  trench. 

The  latter  case  presents  the  most  severe  attack,  albeit  the  least 
likely.  Airblast  loading  could  be  particularly  acute  if  the  trench 
acts  as  a  shock  tube  thereby  extending  (in  the  absence  of  signifi¬ 
cant  attenuation)  the  range  to  a  given  overpressure  beyond  the 
usual  spherical  shock  range.  Thus,  it  is  important  to  under¬ 
stand  the  propagation  of  a  blast  wave  in  the  trench  to  insure  that 
airblast  does  not  impost  excessive  loads  on  the  missile. 

A  number  of  phenomena  will  attenuate  the  blast  wave  as' it 
propagates  in  the  trench.  Among  these  are  wall  ablation,  heat 
transfer  and  momentum  (shear)  loss  to  the  walls,  shock  inter¬ 
actions  with  trench  wall  ribs,  wall  expansion  and  venting. 

Trench  airblast  environments  were  predicted  by  numerical  code 
calculations  primarily  by  Allen  Kuhl  of  TRW  (Reference  1)  and 
Charles  Needham  of  AFWL  (Reference  2).  Kuhl  modeled  the  effects 
of  wall  shear,  expansion  and  venting  in  a  pseudo  one-dimensional 
calculation  while  Needham  performed  two-dimensional  calculations 
of  the  shock  interactions  with  internal  ribs. 

Models  and  assumptions  in  the  codes  must  be  verified  by 
experiment  to  gain  confidence  in  the  predictions.  One-dimensional 
calculations  are  convenient  and  inexpensive  to  run  but  ignore 
the  physics  involved  in  two-dimensional  shock  interactions  and 
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thin  boundary  layers.  Two-dimensional  finite  difference  calcula¬ 
tions  on  the  other  hand  do  not  include  the  effects  of  viscosity. 
Since  both  approaches  were  being  used  to  predict  the  airblast 
environment  in  the  trench,  an  experimental  program  was  undertaken 
as  a  baseline  for  evaluating  the  two  methods. 

Experimental  data  for  shock  propagation  in  ribbed  tubes 
was  nonexistent  in  the  literature  at  the  outset  of  this  study, 
and  only  limited  data  were  available  for  sandgrain  roughened 
tubes  at  lower  shock  strengths  than  of  interest  to  MX.  In  the 
formulation  phase  of  the  present  program,  some  square  rib  shock 
attenuation  tests  performed  by  TRW  and  Martin  Marietta  (Reference 
3)  came  to  our  attention.  Questions  as  to  the  adequacy  of  these 
experiments  were  addressed  in  the  present  study. 

1 . 1  Objectives 

The  present  program  was  undertaken  to  provide  shock  atten¬ 
uation  test  data  in  ribbed  tubes  for  evaluation  of  computer  code 
predictions.  Specific  tests  were  designed  to  investigate  the  in¬ 
fluence  of  Reynolds  number,  shock  strength  and  rib  design  on  shock 
attenuation.  End  wall  reflected  pressure  measurements  were  made 
to  determine  the  effect  of  the  ribs  on  the  blast  plug  loads. 

As  the  program  progressed,  it  became  apparent  that  addi- 
tional  diagnostic  measurements  were  necessary  to  resolve  issues 
surfacing  in  the  airblast  modeling  community.  Consequently, 
the  program  was  expanded  to  measure  the  wall  drag  and  to  provide 
an  optical  description  of  the  flowfield  behind  the  shock. 

1.2  Scope 

A  shock  tube  test  program  was  conducted  at  the  NASA  Ames 
Research  Center  Electric  Arc  Shock  Tube.  To  aid  in  test  design, 
a  one-dimensional  shock  propagation  computer  code  was  developed. 
Test  sections  were  fabricated,  instrumentation  (including  custom 
built  diagnostics  for  wall  drag  and  optical  measurements)  was 
assembled,  and  approximately  65  test  runs  were  conducted. 
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2.0  DESIGN  REQUIREMENTS 

M-X  trench  vulnerability  to  nuclear  attack  was  dictated 
by  the  need  to  provide  adequate  protection  for  the  missile 
against  the  airblast  environment.  Vulnerability  assessments 
require  the  confident  prediction  of  airblast  propagation  within 
the  trench,  and  a  test  program  was  undertaken  to  provide  a 
baseline  for  evaluation  of  computer  code  models. 

2 . 1  M-X  Trench  Configuration 

The  proposed  M-X  trench  system  consists  of  a  number  of 
unmanned  mobile  transporters  making  random  concealed  movements 
in  long  trenches  (tubes)  below  ground.  Airblast  protection 
within  the  trench  is  achieved  by  means  of  blast  plugs  carried 
along  by  the  transporter,  and  the  design  of  these  plugs  is  a 
pacing  item  for  the  trench  concept.  The  baseline  M-X  trench  design 
consists  of  a  4  meter  (13  feet)  diameter  fiber  reinforced  con¬ 
crete  cylinder  with  0.15  meter  (0.5  feet)  high  internal  ribs 
spaced  1.5  meters  (5  feet)  center- to-center  (Figure  1).  The 
ribs  provide  reaction  support  for  the  airblast  plug  as  well  as 
attenuation  of  the  airblast.  Missile  launch  is  achieved  by 
erecting  the  missile  canister  through  the  roof,  and  the  concrete 
cylinder  has  longitudinal  and  transverse  joints  to  ease  the 
breakout  process.  The  joints  provide  strength  only  in  com¬ 
pression;  the  roof  is  held  in  place  by  the  mass  of  the  concrete 
wall  and  the  soil  overburden. 

2 . 2  Shock  Propagation  Analyses 

System  performance  studies  indicate  that  the  blast  plug  must 
survive  a  shock  environment  of  4  x 10®  Pascals  (600  psi)  which  occurs 
for  a  1  MT  surface  burst  at  a  range  of  565m  (1800  feet).  In  the 
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strong  shock  region  near  a  nuclear  burst  on  or  in  the  trench 
massive  wall  ablation  produces  significant  shock  attenuation. 
Wall  expansion  and  venting  are  important  contributors  to  shock 
attenuation  at  high  and  intermediate  shock  strengths.  At  lower 
shock  strengths  (say,  Ap  -  6.8  x  10®  Pascals),  shock  inter¬ 
actions  with  the  ribs  further  reduce  the  shock  strength.  High 
shock  strength  attenuation  is  being  addressed  in  other  DNA 
analytical  and  experimental  programs.  The  present  discussion  is 
concerned  with  the  latter  case,  (i.e.  ,  shock  propagation  in  a 
rigid  tube  with  internal  ribs) . 

2.2.1  Empirical  Prediction  Methods 

A  number  of  shock  tube  experiments  had  been  conducted  at 
low  shock  strengths  in  sand -roughened  tubes  primarily  to  in¬ 
vestigate  attenuation  mechanisms  for  air  entrainment  systems. 
Simplified  models  were  developed  and  used  to  correlate  the  test 
data,  perhaps  the  most  well  known  being  due  to  Porzel  (Reference 
4) .  Skjeltorp  offered  a  more  elaborate  empirical  formulation 
(Reference  5) ,  and  Schiffman  (Reference  3)  correlated  the  TRW 
data  with  a  modified  version  of  these  approaches. 

The  problem  with  these  empirical  models  (as  with  most 
such  approaches)  is  that  the  models  were  fine  tuned  to  match  a 
specific  set  of  data,  and  are  unsuccessful  when  addressing 
problems  outside  the  range  of  parameters  for  which  they  were 
formulated.  Rough  wall  attenuation  of  shocks  in  a  tube  is  far 
from  a  trivial  problem.  Earlier  modeling  efforts  were  too 
simplified  and  led  Kriebel  (Reference  6)  to  conclude,  for  ex¬ 
ample,  " -  from  the  foregoing  comparison  of  theories  -  the 

predicted  attenuation  of  a  shock  front  due  to  tube  wall  friction 
differs  widely  for  both  smooth-walled  tubes  and  rough-walled 
tubes . " 


2.2.2  One-Dimensional  Computational  Models 

Shock  propagation  models  have  been  assembled  by  Kuhl 
at  TRW  (Reference  1)  and  as  part  of  this  program  which  assume 
that  the  flow  is  one-dimensional  and  quasi-steady  (fully  developed 
pipe  flow) .  The  shock  attenuation  mechanisms  invoked  are  wall 
friction  and  heat  transfer.  Solutions  are  obtained  by  numerically 
solving  the  mass,  momentum  and  energy  conservation  equations. 

The  SAI  one-dimensional  model  is  embodied  in  the  RIST 
code  (Reference  7)  which  was  used  to  make  pretest  predictions. 

The  code  solves  the  appropriate  conservation  equations  using 
the  FLIC  method,  and  friction  and  convective  heat  transfer  laws 
are  obtained  from  steady  flow  models.  Schlichting  (Reference  8) 
presents  the  friction  coefficients  (X  =  4Cf)  of  Nikuradse  (Fig¬ 
ure  2)  which  indicate  that  for  high  pipe  Reynolds  number  ReD  and 
large  sandgrain  roughness  kg,  the  friction  coefficient  C^'is 
independent  of  Reynolds  number.  Under  MX  prototype  conditions 
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(ReD  'V/  10  and  R/k  %  13)  the  wall  resistance  formula  for  the 
completely  rough  regime 


(1) 


would  be  expected  to  apply.  For  the  transverse  ribs  in  the  MX 
trench,  an  equivalent  sandgrain  roughness 


(2) 


and  R  is  determined  from  empirical  curves  of  Dalle-Donne  and 
Meyer  (Reference  9) . 
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Figure  2.  Resistance  formula 


Heat  transfer  to  the  wall  is  defined  in  terms  of  a  Stanton 
number  for  which  Dirling's  compressible  flow  model  was  chosen 
(Reference  10) .  Expressed  in  terms  of  a  modified  Reynolds  ana¬ 
logy 

St 

ic“/2) 

where 

+  WTw/pw 

h+  =  -s  —  w  (4) 

ws  v 

w 

is  a  roughness  Reynolds  number  evaluated  at  the  wall. 

Whether  a  one-dimensional  model  is  adequate  depends  to 
a  great  extent  on  whether  the  wall  boundary  closes  the  tube  (fully 
developed  pipe  flow) .  According  to  Mirels  (Reference  11)  ,  the 
shock  boundary  layer  for  a  smooth  wall  MX  trench  requires  a  dis¬ 
tance  of  140  diameters  whereas  a  smooth  wall  shock  tube  of  0.1 
meter  diameter  requires  50  diameters.  Roughness  will  increase 
the  boundary  layer  growth  rate,  but  the  effect  of  large  protuber¬ 
ances  such  as  the  MX  ribs  on  the  boundary  layer  is  unknown. 

2.2.3  Two-Dimensional  Computational  Models 

The  two  dimensional  finite  difference  HULL  code  (Ref¬ 
erence  12)  has  been  used  by  Needham  at  the  Air  Force  Weapons 
Laboratory  to  calculate  the  airblast  environment  in  an  MX  trench. 
HULL  solves  the  finite  difference  analogs  of  a  system  of  partial 
differential  equations  describing  an  inviscid,  nonconducting  fluid. 
A  two-dimensional  cylindrical  geometry  was  used  to  represent  the 

7 

MX  trench,  and  a  6.8  x  10  Pascal  overpressure  shock  waveform  from 
a  1  MT  surface  burst  was  input  at  the  trench  entrance  (Reference  2) . 

The  HULL  code  calculation  revealed  a  complicated  two 
dimensional  flow  structure  due  to  shock  interaction  with  the 
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ribs,  and  the  shock  decayed  more  rapidly  than  an  equiv¬ 
alent  free  air  burst.  Shock  attenuation  was  caused  by  the 
nonisentropic  conversion  of  kinetic  energy  to  thermal  energy  and 
the  redistribution  of  axial  momentum  into  radial  components.  The 
influence  of  viscosity  on  these  calculations  has  not  been  deter¬ 
mined. 


2 . 3  Experimental  Requirements 

Since  the  one-dimensional  and  two-dimensional  models  treat 
different  physical  mechanisms,  it  was  necessary  to  provide  a  means 
of  evaluating  the  approaches.  The  philosophy  adopted  was  to 
obtain  a  set  of  experimental  data  for  a  well  described  flowfield 
and  scaled  geometry  as  a  test  for  the  computer  codes.  Thus,  if 
a  code  could  adequately  predict  the  experimental  results,  confi¬ 
dence  would  be  gained  in  the  prototype  predictions. 

Small  scale  experiments  can  be  directly  useful  to  proto¬ 
type  assessments  if  it  can  be  shown  that  the  scale  either  does 
not  influence  the  test  results  or  can  be  accounted  for.  For  a 
geometrically  scaled  test  model,  the  model  size  enters  through 
the  nondimen sional  Reynolds  number.  Skin  friction,  for  example, 
is  usually  Reynolds  number  dependent;  however,  for  rough  tubes 
at  sufficiently  high  Reynolds  number,  the  skin  friction  is 
Reynolds  number  independent.  Thus,  determining  the  influence 
(or  lack  thereof)  of  Reynolds  number  was  a  test  requirement. 

Shock  attenuation  rates  were  found  to  be  a  function 
of  shock  strength  by  TRW  and  Martin  for  tubes  with  internal  ribs 
and  by  Porzel  for  sandgrain  roughened  tubes.  Therefore,  a  re¬ 
quirement  of  the  experiments  was  to  perform  the  tests  from  as 
high  an  initial  shock  strength  as  possible  down  to  an  initial 
shock  strength  of  40  (equivalent  of  a  600  psi  shock) .  To  main¬ 
tain  high  Reynolds  numbers  and  absolute  shock  pressures,  tests 
were  to  be  conducted  as  near  ambient  initial  pressure  as  prac¬ 
tical.  Since  the  ribs  were  expected  to  perturb  the  flow,  end 
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wall  reflected  pressure  measurements  were  required  to  estimate 
the  load  history  on  the  M-X  blast  plug. 

The  primary  rib  geometry  for  the  experiments  was  to 
be  the  nominal  MX  design  (Figure  1).  However,  it  was  also 
of  interest  to  determine  if  the  rib  configuration  could  be  opti¬ 
mized  to  provide  increased  attenuation.  The  empirical  skin  fric¬ 
tion  correlation  of  Dalle  Donne  and  Meyer  indicated  that  the  M-X 
rib  separation  distance  was  well  chosen  but  that  increased  drag 
could  be  achieved  (in  steady  flow)  by  increasing  the  height  to 
width  ratio.  A  second  rib  geometry  was  built  to  test  this 
hypothesis. 

Details  of  the  shock  wave  flowfield  were  deemed  necessary 
to  understand  the  physics  of  the  shock  wave  attenuation  process. 

For  the  one-dimensional  models,  the  question  of  appropriate  wall 
drag  models  and  the  state  of  the  boundary  layer  were  controversial 
and  required  experimental  investigation.  In  addition,  a  qualita¬ 
tive  picture  of  the  rib/f lowf ield  interaction  was  necessary  to  guide 
the  development  of  one  and  two-dimensional  models. 
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3.0 


SHOCK  TUBE  TEST  PROGRAM 


A  series  of  shock  tube  experiments  were  conducted  at  the 
NASA  Ames  Research  Center  Electric  Arc  Shock  Tube  (EAST)  to  aid 
in  the  understanding  of  airblast  propagation  in  the  MX  trench. 

The  facility  was  selected  because  of  its  ability  to  produce  strong, 
repeatable  shock  waves  at  ambient  or  near  ambient  initial  pres¬ 
sures  . 

Test  sections  containing  removable  liners  with  internal 
ribs  were  attached  to  the  EAST  facility  and  tested  at  three  initial 
shock  strengths  after  facility  calibration  tests.  Blast  plug  load 
histories  were  obtained  from  pressure  measurements  on  a  rigid  end 
wall  located  at  the  end  of  a  ribbed  test  section.  The  nominal 
instrumentation  consisted  of  time  of  arrival  (TOA)  probes  and  high 
frequency  pressure  transducers.  In  addition,  custom  built  diagnos¬ 
tics  were  developed  to  determine  the  loads  on  a  floating  wall  ele¬ 
ment  and  to  provide  optical  measurements  of  the  rib/f lowfield  inter¬ 
actions. 


3 . 1  Facility  Description 

The  NASA  Ames  Research  Center  EAST  facility  (Figure  3)  is 
an  electric  arc  shock  tube  used  for  high  energy  shock  physics  ex¬ 
periments.  EAST  consists  of  a  1.2  MJ  capacitor  bank,  an  0.1  meter 
I.D.  shock  tube  and  a  dump  tank  (Figure  4).  For  the  present  tests 
a  1.36  meter  long  cylindrical  driver  was  employed. 

Strong  shocks  are  created  by  discharging  energies  approach¬ 
ing  1  MJ  into  a  driver  section  through  an  exploding  wire  techni¬ 
que.  Joule  heating  of  the  driver  gas  takes  place  when  an  arc  is 
struck  along  the  wire  path  between  the  ground  electrode  at  the 
diaphragm  station  and  the  high  voltage  electrode  at  the  base  of 
the  driver  (Figure  5) .  Alternate  gases  such  as  helium  are  used 
in  the  driver  section  to  create  stronger  shocks  than  would  occur 
with  air  in  the  driver.  Details  of  the  driver  operation  can  be 
found  in  References  13  and  14. 
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Figure  4.  NASA  Ames  electric  arc  shock  tube  schematic. 
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Operation  of  the  shock  tube  is  easily  illustrated  by  con¬ 
sidering  a  constant  diameter  tube  divided  into  two  regions:  a 
high  pressure  driver  section  and  a  lower  pressure  driven  section 
separated  by  a  diaphragm  (Figure  6a) .  When  the  diaphragm  is 
burst,  a  shock  wave  propagates  down  the  driven  section,  and  an 
expansion  fan  propagates  back  into  the  driver  section  (Figure 
6b) .  The  expansion  fan  converts  the  stagnation  energy  of  the 
driver  gas  into  flow  energy  of  the  expanded  driver  gas  and  the 
work  done  in  compressing  and  accelerating  the  driven  gas.  Ex¬ 
pansion  waves  travel  at  the  speed  of  sound  relative  to  the  local 
driver  gas  conditions;  thus  the  driver  gas  sound  speed  is  an 
improtant  paramerter  in  creating  strong  shocks.  References  15 
and  16  present  the  details  of  shock  tube  performance. 

Using  the  notation  of  Figure  6,  the  ideal  shock  Mach  num¬ 
ber  is  determined  from 
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where  A  is  sound  speed  and  the  double  subscript  implies  a  ratio. 
The  shock  pressure  ratio  or  strength  is 
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Driver  gas  (Region  4)  state  after  arc  discharge  is  estimated  from 
constant  volume  heating  for  ideal  gases 
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where  m  is  the  molecular  weight.  Equations  5  through  8  were  used 
to  estimate  the  shock  tube  performance;  final  test  conditions 
(Figure  7)  were  determined  from  calibration  runs. 

Velocity  and  state  variables  behind  the  shock  wave  (Region 
two)  are  found  from  the  initial  state  of  the  driven  tube  gas  and 
the  shock  Mach  number 
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Reynolds  number  based  on  diameter  in  region  two  is 
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where  STP  refers  to  standard  conditions,  and  the  viscosity  is  taken 
as  proportional  to  temperature. 

Since  velocity  and  pressure  are  constant  across  the  contact 
surface,  conditions  in  Region  three  are  found  from  the  initial  state 
of  the  driver  gas  and  the  velocity  in  Region  3. 
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The  Region  3  Reynolds  number  is 
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Equations  12  and  16  were  used  to  estimate  the  flow  Reynolds  number 
for  scaling  considerations. 

Since  the  contact  surface  follows  closely  behind  the  shock 
wave  for  high  shock  strength  flows,  Region  3  gases  occupy  the  test 
section  for  a  large  portion  of  the  test.  Communication  of  the 
effects  of  the  ribs  on  Region  3  gases  to  the  shock  is  con¬ 
trolled  by  the  acoustic  impedance  ratio  across  the  contact  sur¬ 
face 


r32  p 32  A32 
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If  T  is  near  one  acoustic  signals  are  transmitted  intact,  if  T 
is  largely  different  from  one  the  signals  can  be  weakened  (T  >  1) 
or  amplified  (T  <  1) . 

3 . 2  Models  and  Instrumentation 

Two  test  sections  were  constructed  for  airblast  attenuation 
experiments.  Cylindrical  4142  steel  seemless  tubes  0.18  meter 
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(7  inch)  O.D.  by  0.13  meter  (5  inch)  I.D.  were  bolted  to  the  EAST 
facility  via  existing  flanges  threaded  onto  the  steel  tubing. 

Figure  8  shows  a  single  3.1  meter  (10  feet)  long  test  section  in 
place  for  the  initial  attenuation  experiments;  the  final  experi¬ 
ments  were  run  with  a  5.5  meter  (18  feet)  test  section. 

Interchangable  liners  with  internal  ribs  were  inserted 
into  the  shock  tube  test  sections.  The  internal  diameter  of  the 
liners  (rib  tops)  matched  the  NASA  shock  tube  smooth  wall  diameter 
(10  cm)  to  avoid  an  abrupt  area  change.  Two  sets  of  liners  were 
tested  (Figure  9) :  one  set  was  a  1/39  scale  (based  on  diameter) 
replica  of  the  proposed  MX  rib  geometry  (referred  to  as  MX) ,  and 
one  set  was  an  attempt  to  optimize  shock  attenuation  via  the 
Dalle  Donne  and  Meyer  skin  friction  correlations  within  the  MX 
rib  height  constraints  (referred  to  as  DD&M) . 

Pressure  measurements  were  made  along  the  test  section  pri¬ 
marily  at  the  rib  tops;  a  few  rib  valley  measurements  indicated  no 
significant  differences.  Pressure  transducers  were  mounted  in 
steel  inserts  in  the  sidewall  to  protect  them  from  liner/wall 
relative  motion.  Because  the  mountings  were  wider  than  the  in¬ 
dividual  ribs,  the  space  between  two  ribs  was  filled  in  at  the 
trandsucer  locations  (Figure  10) . 

Several  initial  tests  were  conducted  to  provide  informa¬ 
tion  for  an  MX  interchange  meeting;  the  test  arrangement  con¬ 
sisted  of  2.4  meters  (8  feet)  of  MX  ribs  followed  by  0.6  meters 
(2  feet)  of  smooth  wall.  The  MX  ribs  began  at  the  shock  tube 
interface  (Figure  11a),  and  the  first  pressure  station  was  0.3 
meters  (1  foot)  from  the  first  rib.  The  initial  tests  indicated 
a  dramatic  drop  in  shock  strength  by  the  time  the  shock  had 
traversed  the  8  ribs  to  the  first  pressure  transducer.  To 
investigate  this  entrance  region  effect  the  second  test  section 
was  modified,  and  the  remainder  of  the  tests  were  conducted  with 
a  pressure  transducer  downstream  of  the  first  rib  (Figure  lib) . 
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Figure  10,  Pressure  transducer  mounting. 


Figure  11a.  Initial  entrance  region  design. 


Roughness  Begins  (X  =  4.0  meters) 


Station  R1 


Figure  11b.  Final  entrance  region  design. 
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Shock  time  of  arrival  (TOA)  measurements  were  made 
along  the  shock  tube  as  a  basis  for  code  comparison  and  as 
a  check  on  the  pressure  measurements.  TOA  measurements  were 
accomplished  with  NASA  twin  probe  ionization  gauges  coupled  to 
a  counter  system  via  existing  electronics.  The  probes  operate 
on  the  principle  that  the  heated  gas  behind  the  shock  front  pro¬ 
vides  a  conductive  path  between  the  probes  which  produces  a  trig¬ 
ger  signal.  TOA  locations  are  given  in  Appendix  A. 

Pressure  history  measurements  were  made  at  various 
stations  along  the  smooth  and  ribbed  test  sections  (Figure  12) . 
Smooth  wall  stations  were  designated  C  through  I;  ribbed  wall 
stations  were  designated  Rl  through  R6  for  transducers  mounted 
flush  with  the  rib  tops  and  V2  through  V3  for  transducers  mounted 
in  the  rib  valleys.  Transducer  locations  for  the  final  runs 
are  given  in  Appendix  B.  For  the  blast  plug  tests,  the  pressure 
transducers  were  mounted  at  three  radial  stations  (r  =  0,  R/2 
and  3R/4  in  a  rigid  end  wall  located  downstream  of  a  24  diameter 
long  ribbed  section.  PCB  Piezotronics  quartz  pressure  trans¬ 
ducers  models  113A22  and  113A24  provided  high  frequency  (1  ysec 
rise  time)  measurements.  The  transducers  were  protected  from  the 
thermal  environment  by  covering  the  transducer  face  with  a  thin 
layer  of  RTV. 

A  custom  built  load  cell  (PCB  Model  217M03)  was  employed 
to  measure  the  longitudinal  loads  on  a  0.15  meter  (0.5  feet)  long 
floating  M-x  rib  liner  section  containing  four  ribs.  Liner  sec¬ 
tions  fore  and  aft  of  the  floating  liner  were  pinned  in  place,  and 
a  nylon  ring  was  loaded  in  compression  in  front  of  the  liner  to 
take  up  any  relative  motion  (Figure  13) .  The  load  cell  had  to 
be  specially  designed  and  fabricated  since  none  could  be  found 
with  an  0.1  meter  (4  inch)  I.D.  to  allow  the  shock  to  propagate 
through  the  test  section.  Fast  response  was  required  to  measure 
the  shock  induced  loads.  Based  on  a  measured  stiffness  of  3  x  lO10 
newton/meter  and  a  liner  mass  of  4.9  kilograms,  the  load  cell/liner 
response  time  was  estimated  to  be  of  the  order  of  80  usee  for  a 
simple  spring-mass  system. 
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A  square  test  section  was  constructed  to  accomodate 
optical  measurements  of  two-dimensional  rib/f lowf ield  inter¬ 
actions  (Figure  14) .  To  avoid  complicated  geometrical  transi¬ 
tions  from  the  cylindrical  NASA  shock  tube,  a  cookie  cutter 
approach  was  used,  and  this  limited  the  optical  test  section 
internal  dimension  to  5  cm  (2  inches).  Scaled  MX  ribs  (scale 
factor  78)  were  located  on  the  top  and  bottom  of  the  test  sec¬ 
tion  extending  67  cm  (26  inches)  upstream  and  30  cm  (12  inches) 
downstream  of  the  center  of  a  5  cm  x  15  cm  optical  port  utiliz¬ 
ing  plexiglass  windows  flush  mounted  with  the  sidewalls  (Figure 
15).  The  overall  test  section  length  was  151  cm  (5.0  feet) 
extending  30  cm  (12  inches)  into  the  forward  NASA  shock  tube  and 
ending  abruptly  at  the  aft  NASA  shock  tube.  Pressure  trans¬ 
ducers  were  located  at  both  ends  of  the  test  section  (labeled 
TS1  and  TS3)  and  in  the  optical  port  region  (TS2) . 

An  SAI  custom  built  holographic  interferometer  was 
employed  to  record  shadowgrams  and  interferograms  of  the  flow- 
field.  A  Nd/YAG  frequency  doubled  laser  emitting  30  mjoules  of 

o 

5320  A  light  in  20  nsec  pulses  provided  a  source  of  coherent 
illumination.  Agfa-Gavaet  film  (number  10E-56)  which  spatially 
resolves  1500  lines/mm  was  used  to  record  the  images.  Figure  16 
shows  the  optical  arrangement  in  place  at  the  NASA  shock  tube; 
details  of  the  holography  system  are  contained  in  Appendix  C. 

An  ionization  gauge  at  Station  I  detected  the  shock 
arrival,  and  the  gauge  output  triggered  a  delay  pulse  generator 
which  fired  the  laser  after  a  selected  time.  The  actual  laser 
firing  time  was  verified  by  a  detector  placed  within  the  laser 
beam.  Shock  wave  arrival  along  the  test  section  was  measured  at 
three  pressure  transducer  stations  relative  to  shock  wave  TOA  at 
Station  I.  Therefore,  the  time  at  which  the  laser  fired  could  be 
determined  relative  to  the  shock  wave  arrivel  in  the  optical  port 
for  comparison  with  code  calculations  and  to  estimate  the  posi¬ 
tion  of  the  shock  relative  to  the  optical  port. 
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3.3 


Data  Recording  and  Processing 

Shock  wave  TOA  measurements  were  made  in  the  smooth 
sections  using  the  NASA  MSTIMS  (Multi-Station  Time  Interval 
Measuring  System)  which  is  a  ten  channel  counter  system  with  0.1 
microsecond  time  resolution  (Reference  17).  Counting  begins 
after  the  arc  discharge  is  extinguished  (without  this  precaution 
the  counters  could  be  started  prematurely  by  RF  noise  from  the 
discharge) .  The  counters  are  stopped  by  trigger  signals  from 
the  NASA  twin  probe  ionization  gauges  located  along  the  shock 
tube.  The  outputs  of  the  ionization  gauges  are  monitored  by 
Schmidt  triggers  which  are  enabled  at  an  adjustable  time  after 
arc  discharge  (usually  slightly  longer  than  the  time  delay  to 
begin  counting) .  When  the  Schmidt  trigger  is  fired,  a  50  V  pulse 
is  generated  which  stops  the  counters.  In  the  ribbed  wall  test 
sections  TOAS  were  obtained  from  the  pressure  transducers. 

Pressure  histories  at  various  stations  were  obtained  from 
Polaroid  records  of  the  pressure  transducer  signals  displayed  on 
Tektronix  Model  555  Oscilloscopes.  Vendor  supplied  calibrations 
(Appendix  B)  were  used  to  convert  voltage  to  pressure  levels.  It 
was  estimated  that  the  shock  wave  pressure  ratios  could  be 
determined  from  the  photographs  to  with  ±  10  percent;  shock  wave 
TOA  was  measured  to  ±  5  usee.  Load  cell  data  were  recorded  and 
processed  in  a  similar  manner.  Details  of  interferogram  data 
analysis  are  presented  in  Appendix  C. 

3 . 4  Test  Matrix 

Shock  tube  driver  conditions  were  selected  to  produce 
shock  strengths  ?2i  °f  80,300  and  700  at  the  entrance  to  the 
ribbed  test  section  for  the  shock  attenuation  tests.  Since  the 
test  flowfield  is  largely  composed  of  Region  3  gas,  the  tests 

were  designed  to  keep  the  acoustic  impedance  ratio  r23  across 

the  contact  surface  as  close  to  one  as  practical  to  insure 
that  wall  effects  would  be  communicated  to  the  shock  wave. 
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Table  1  gives  the  conditions  for  the  shock  attenuation 
tests.  Shock  tube  repeatability  as  measured  by  run  to  run  varia¬ 
tion  of  shock  Cg  between  stations  at  2.0  and  3.2  meters  from  the 
diaphragm  was  generally  better  than  5  percent.  Four  series  of 
tests  are  indicated.  The  first  two  series  compared  Reynolds  num¬ 
ber  effects  at  a  constant  shock  strength  (P2i  =  80).  Series  three 
and  four  addressed  the  effect  of  increased  shock  strength  (P2i  = 
300  and  700)  at  nearly  constant  Reynolds  number.  Runs  5,  10, 

7  and  50  were  smooth  tube  calibration  tests.  Runs  56,  54  and  52 
enployed  DD&M  ribs  while  the  remainder  tested  MX  ribs.  Runs  58 
and  59  demonstrated  test  repeatability.  End  wall  reflected  pres¬ 
sures  were  measured  on  Run  63. 

Test  conditions  for  the  optical  runs  were  the  same  as  for 

Series  2,  except  that  the  driver  pressure  was  raised  to  5.4 

atmospheres.  The  shock  wave  pressure  ratio  in  the  optical  port 

(Station  TS2 )  was  P22  =  70  (Appendix  B) ,  and  the  mean  shock 

speed  between  stations  TSl  and  TS2  was  2.60  Km/sec  (M  =  7.47) 

s 

with  a  variation  of  two  percent  for  Runs  37  through  49.  The 
optical  system  configuration  for  each  run  and  firing  time  re¬ 
lative  to  the  station  TS2  TOA  are  listed  in  Table  2. 
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Table  1.  Shock  attenuation  test  conditions. 


Shock  velocity  between  stations  D  and 


e  2.  Optical  test  series. 


4 . 0  EXPERIMENT  RESULTS 

Sixty  five  runs  were  made  in  the  NASA  East  facility  to 
investigate  the  propagation  of  strong  shock  waves  in  tubes  with 
internal  ribs.  Calibration  runs  established  the  shock  wave  char¬ 
acteristics  in  the  existing  smooth  wall  NASA  shock  tube  and  deter¬ 
mined  the  arc  driver  performance.  A  few  attenuation  tests  were 
conducted  using  one  ribbed  test  section  with  an  L/D  of  24  to  pro¬ 
vide  initial  data  in  time  for  an  MX  interchange  meeting.  Final 
attenuation  data  and  wall  drag  measurements  were  obtained  in 
two  test  sections  with  a  combined  L/D  of  48.  Only  these  latter 
tests  are  reported  here.  End  wall  reflected  pressures  were 
obtained  downstream  of  the  first  test  section  and  optical  data  were 
obtained  in  the  two  dimensional  test  section. 

Test  results  are  presented  in  their  reduced  form  in  the 
following  section.  The  raw  TOA  data  and  pressure  traces  are  cata¬ 
logued  in  Appendixes  A  and  B  along  with  appropriate  transducer 
sensitivities  and  oscilloscope  settings. 

4 . 1  Pressure  Waveforms 

A  conventional  shock  tube  produces  a  shock  wave  which  is 
followed  by  a  region  of  constant  pressure.  The  constant  pressure 
region  is  terminated  by  the  front  of  the  rarefaction  wave  reflecting 
from  the  aft  end  of  the  driver.  The  rate  at  which  the  rarefaction 
front  approaches  the  shock  and  at  which  the  shock  decays  after  the 
front  reaches  the  shock  are  dependent  on  the  sound  speed  of  the 
driver  gas  and  the  length  of  the  driver.  Thus,  pressure  waveforms 
can  be  tailored  to  a  certain  extent  by  controlling  these  two  vari¬ 
ables  . 

For  the  present  tests,  it  was  desired  that  the  shock  wave 
traverse  the  entire  length  of  the  test  section  (4.8  meters)  before 
rarefaction  waves  reached  the  entrance  in  order  to  investigate  rib 
effects  in  the  simplest  possible  flowfield.  The  longest  available 
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driver  (1.36  meters)  was  used  to  maximize  the  constant  pressure 
regime.  Pressure  waveforms  (e.g.  Figure  17)  in  the  smooth  shock  tube 
calibration  runs  exhibited  a  constant  pressure  duration  which 
depended  on  shock  strength  and  distance  along  the  shock  tube.  In 
order  of  increasing  shock  strength  the  constant  pressure  durations 
were  100,  600,  and  900  ysec  at  the  entrance  to  the  ribbed  test 
section  (L/D  =  40) . 


Figure  17,  Sample  pressure  waveform  (26.2  atm/div 
and  500  j/sec/div)  , 


The  effect  of  MX  ribs  on  the  pressure  waveform  for  P2^  =  80 
was  to  stretch  out  the  constant  pressure  region  and  reduce  the 
pressure  ratio  as  the  shock  traveled  through  the  ribbed  test  sec¬ 
tion  (Figure  18) .  However,  when  the  P2^  =  300  shock  encountered 
the  ribbed  test  section  the  pressure  waveforms  exhibited  a  gradual 
rise  behind  the  shock  wave  to  two  or  three  times  the  local  shock 
value  (Figure  18) .  Pressure  waveforms  for  the  strongest  incident 
shocks  tested  (P21  =  700)  showed  a  similar  behavior.  The  increasing 
pressure  may  not  be  of  particular  significance  to  MX  trench  design, 
as  trench  expansion  and  venting  effects  (not  simulated  in  these 
tests)  will  tend  to  reduce  the  pressure  rise.  However,  calcula¬ 
tions  of  the  experiments  can  only  be  considered  successful  if  these 
features  of  the  pressure  waveforms  are  predicted. 
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Pressure  waveforms  in  the  smooth  tube  forward  of  the 


ribbed  test  section  showed  several  discontinuities  well  after  the 
rarefaction  waves  had  entered  the  test  section  indicating  shocks 
were  moving  upstream  into  the  smooth  shock  tube  section.  In  the 
TRW  (Reference  3)  tests  these  perturbations  were  thought  to  be 
due  to  reflections  off  the  ribs  projecting  into  the  flow;  how¬ 
ever,  in  the  present  tests  the  rib  tops  were  aligned  with  the 
smooth  tube.  One-dimensional  calculations  indicated  that  the 
compression  region  forming  in  the  ribbed  tube  builds  into  a  rela¬ 
tively  weak  shock  wave  which  propagates  back  upstream  as  the  up¬ 
stream  pressure  decays.  Two  dimensional  calculations  are  required 
to  resolve  the  detailed  flow  processes. 

The  waveforms  presented  here  were  all  measured  with  trans¬ 
ducers  located  in  the  rib  tops.  No  significant  differences  were 
noted  in  the  few  valley  measurements  (e.g..  Figures  B-7  and 
B-8  of  Appendix  B) . 

4 . 2  Shock  Attenuation 

From  the  x-t  plots  prepared  from  the  TOA  data  (Figure  19 
to  21) ,  it  is  apparent  that  the  shock  attenuates  in  the  ribbed 
test  sections.  The  P21  =  80  shock  traveled  48  diameters  (4.8 
meters)  in  the  smooth  tube  in  1.89  msec.  In  the  MX  ribbed  test 
section  the  traversal  time  was  2.30  msec,  and  in  the  DD&M  test 
section  the  traversal  time  was  2.67  msec.  For  the  P2^  =  700  shock 
the  comparable  times  are  0.65  msec  in  the  smooth  section  and  0.87 
msec  in  the  MX  ribbed  section.  Thus,  the  rarefactions  from  the 
driver  entered  the  test  section  before  the  shock  wave  left  for  the 
P21  =  runs.  For  the  P21  =  700  runs  the  shock  wave  left  the 
test  section  before  the  rarefactions  arrived  at  the  entrance. 

Also  apparent  in  the  x-t  diagrams  are  the  trajectories 
of  the  shocks  propagating  upstream.  Assuming  that  there  is  little 
or  no  flow  down  the  shock  tube  (the  incoming  pressure  waves 
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Figure  20.  Shock  TOA  for  DD&M  ribs  (Poi  =  80). 


TIME  (nsec) 


Figure  21.  Shock  TOA  for  MX  ribs  (Pot  =  700 


have  decayed  due  to  rarefactions  from  the  driver) ,  the  shallow 
trajectories  provide  further  evidence  that  the  returning  shocks 
are  relatively  weak. 

Pressure  measurements  on  the  rib  tops  also  demonstrate  the 
shock  attenuation.  For  =  80  with  an  atmospheric  test  section, 

the  shock  attenuates  to  about  76  percent  of  its  smooth  wall  value 
by  24  diameters  of  travel  and  to  45  percent  of  its  smooth  wall 
value  by  48  diameters  (Figure  22).  Three  runs  (57,  58,  59)  indi¬ 
cate  the  repeatability  of  the  pressure  measurements;  the  data 
generally  agree  within  ±  10  percent  (Figure  22) .  Shock  pressure 
ratios  from  the  TOA  measurements  via  the  Rankine-Hugoniot  relation¬ 
ships  were  generally  within  15  percent  of  those  determined  from 
pressure  gauges  located  midway  between  the  TOA  probes.  For  ex¬ 
ample,  the  measured  shock  pressure  ratio  at  the  X/D  =  20  station 
in  Figure  22  was  95.3  while  the  value  determined  from  TOA  measure¬ 
ments  was  83.5.  Consequently,  for  the  remainder  of  the  test  series 
only  one  successful  (defined  as  acquiring  the  necessary  data) 
test  run  was  conducted  for  each  objective,  and  only  the  reduced 
pressure  data  are  discussed. 

4.2.1  Reynolds  Number  Effect  on  Shock  Attenuation 

Shock  attenuation  data  taken  from  the  pressure  measure¬ 
ments  for  ■  =  SO  showed  no  discernable  effects  of  flow  Rey¬ 
nolds  number.  Tests  run  with  initial  pressures  of  1.0  and  0.1 
atmospheres  produced  nearly  identical  decay  rates  although  the 
flow  Reynolds  numbers  differed  by  an  order  of  magnitude  (Figure 
23) .  The  absolute  values  of  the  shock  strengths  were  slightly 
different  for  the  two  tests  due  to  sensitivity  of  the  arc  dis¬ 
charge  to  initial  driver  pressure;  however,  a  consistent  trend  is  seen 
in  the  smooth  and  ribbed  wall  data.  In  fact,  if  the  data  is  shifted 
to  match  in  the  smooth  tube,  the  results  overlap. 
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Figure  22.  Shock  attenuation  by  MX  ribs. 
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Figure  23.  Reynolds  number  effect  on  shock  attenuation. 
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The  implication  here  is  that  scale  effects  introduced 
through  the  Reynolds  numbers  are  not  important  to  the  attenua  - 
tion  process.  This  could  come  about  in  two  ways:  the  pipe 
flow  type  skin  friction  is  truly  in  the  Reynolds  number  inde¬ 
pendent  region  or  the  attenuation  mechanism  is  a  purely 
inviscid  shock-rib  interaction  phenomena. 

4.2.2  Shock  Strength  and  Rib  Geometry 

Effects  on  Shock  Attenuation 

Shock  attenuation  data  taken  from  the  pressure  measure¬ 
ments  are  summarized  in  Figures  24  through  26  for  three  initial 
shock  strengths  and  two  rib  geometries.  The  present  data  indicated 
essentially,  a  relatively  uniform,  exponential  decay  rate  with  the 
shock  waves  reaching  half  their  associated  smooth  wall  strength  in 
about  30  diameters  for  =  80  and  in  about  20  diameters  for 

P21  =  and  700- 

It  must  be  noted  that  measurements  made  in  the  preliminary 
portion  of  this  program  indicated  a  dramatic  drop  in  shock  strength 
at  a  station  three  diameters  from  the  rib  test  section  entrance 
(Figure  27) .  The  present  tests  conducted  with  the  first  transducer 
located  behind  the  first  rib  and  the  second  transducer  nine  dia¬ 
meters  from  the  start  of  the  ribs  are  not  conclusive  as  to  how 
quickly  the  shock  decays  initially.  Caution  should  be  used  in 
assuming  simple  exponential  data  fits  until  the  entrance  effect 
is  resolved. 

The  DD&M  rib  geometry  appeared  to  provide  slightly  better 
shock  attenuation  than  the  MX  ribs.  The  increased  attenuation  was 
not  evident  until  the  shock  had  traveled  through  20  to  30  diameters 
of  the  ribbed  section  after  which  the  effect  seemed  to  be  increasing. 
The  additional  attenuation  at  the  end  of  the  test  section  (L/D  =  48) 
was  25,  27  and  19  percent  over  the  MX  attenuation  in  order  of 
increasing  shock  strength. 


SHOCK  STRENGTH  (P0i) 


Figure  24.  Shock  attenuation  at  ?2i  =  80. 
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Figure  26.  Shock  attenuation  at  ?2i  =  700. 
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Figure  27.  Shock  attenuation  from  preliminary  tests 


4.3 


Optical  Flowfield  Measurements 

Holographic  techniques  provided  a  qualitative  description 
of  the  rib/f lowf ield  interaction  as  well  as  limited  quantitative 
measurements  of  the  density  fields.  Photographs  of  the  holo¬ 
graphic  images  are  catalogued  in  Appendix  C;  shadowgrams  and 
interf erograms  of  the  region  near  the  main  shock  and  at  7.5  dia¬ 
meters  behind  the  main  shock  are  described  here.  Additional  quanti¬ 
tative  data  could  be  reduced  from  the  holographic  interferogram 
plates . 

Laser  shadowgrams  of  the  main  shock  region  (Figure  28) 
confirm  the  existence  of  a  multiple  shock  structure  predicted  by 
AFWL  two-dimensional  calculations  (Reference  2) .  The  source  of 
the  multiple  shock  structure  is  the  reflection  of  the  main  shock 
from  the  ribs.  Each  shock-rib  interaction  sends  out  a  blast  wave¬ 
like  shock  which  interacts  with  the  main  shock  and  with  the  other 
secondary  shocks.  This  propagation  controls  the  strength  of  the 
main  shock  and  therefore  drives  the  attenuation  process.  It  is 
useful  to  discuss  schematically  the  shock  rib  interaction  for  one 
rib,  and  then  to  note  the  additional  features  of  multiple  ribs  from 
the  flow  visualization. 

As  the  main  shock  passes  each  trapezoidal  rib  compression 
and  expansion  disturbances  are  created  at  both  the  forward  and 
rearward  faces.  Mach  reflection  occurs  on  the  forward  face 
(Reference  23)  causing  a  reflected  shock  to  propagate  along  the 
backside  of  the  main  shock.  Shock  strength  is  discontinuous 
across  this  three  shock  intersection,  and  the  trajectory  of  this 
point  is  called  a  triple  point  path  (Reference  19).  Mach  reflection 
continues  until  the  Mach  stem  reaches  the  corner  where  an  expansion 
fan  is  formed.  The  expansion  fan  overtakes  and  weakens  the  Mach 
stem  causing  the  curvature  required  to  keep  the  shock  normal  to  the 
top  of  the  rib.  As  the  expansion  fan  overtakes  the  triple  shock 
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igure  28.  Laser  shadowgram  of  shock/rib  interactions. 


the  reflected  wave  is  weakened  reducing  the  shock  strength  dis¬ 
continuity  at  the  triple  shock.  This  weakening  also  decelerates 
the  radial  velocity  of  the  triple  shock  which  in  the  limit  of 
weak  reflected  shocks  is  the  sonic  velocity  behind  the  shock. 

Another  expansion  fan  is  created  when  the  shock  reaches 
the  downstream  edge  of  the  rib.  Again,  the  expansion  fan  over¬ 
takes  and  weakens  the  portion  of  the  main  shock  near  the  wall 
keeping  the  shock  normal  to  the  rear  face  of  the  rib.  Mach  re¬ 
flection  occurs  when  the  shock  on  the  rear  fan  reaches  the  bottom 
of  the  rib  creating  a  second  triple  shock  intersection.  This 
event  is  the  last  portion  of  one  cycle  of  the  motion  of  the  main 
shock. 

The  rib/shock  interactions  combine  in  a  complicated  but 
symmetric  pattern.  The  main  features  observed  in  the  shadowgram 
are  identified  in  Figure  29,  and  they  include: 

•  Reflected  secondary  shocks  on  the  main 
shock  forming  triple  shock  intersections. 

•  Slip  streams  in  the  wake  of  the  triple  shock 
intersections. 

•  Vortex  like  separated  regions  on  the  top  and 
backside  of  the  robs. 

•  Reflected  shock  interactions  with  turbulent 
cavity  flows. 

To  quantify  the  two  dimensional  aspects  of  the  flow, 
holographic  interferometry  was  used  to  measure  density  variations 
across  the  shock  tube  (Appendix  C) .  A  holographic  interferogram 
of  the  main  shock  (Figure  30)  displays  a  complicated  image  of  the 
flow.  The  density  profile  at  section  A-A  one  rib  spacing  behind 
the  shock  (Figure  31)  reveals  a  strong  compression  process  with 
the  density  at  the  wall  exceeding  twice  the  centerline  density. 
This  increased  density  is  presumably  a  result  of  the  reflection 
of  the  main  shock  from  the  ribs. 
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The  main  features  of  the 
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Figure  30,  Holographic  Interferogram  of  shock/rib  interactions. 


SHOCK  (Pi  =  9,03  x  10-4  gm/cm^), 


Farther  behind  the  shock  (X/D  =  2.5)  the  flowfield  appears 
to  become  stationary  with  bow  shocks  extending  diagonally  across 
the  shock  tube  (Figure  32) .  The  angle  between  the  straight  por¬ 
tion  of  the  bow  waves  and  the  shock  tube  axis  is  approximately  the 
Mach  angle  (a  =  30°)  where  the  flow  Mach  number  is  calculated 
from  the  shock  speed  via  the  Rankine-Hugoniot  equations.  This 
would  imply  relatively  weak  waves. 

A  holographic  interferogram  (Figure  33)  of  the  flow  at 
a  later  time  than  the  laser  shadowgram  (X/D  =  7.5)  displays  the 
sensitivity  of  the  technique  to  detect  small  density  changes. 

The  dominant  features  of  the  flow  are  displayed  in  the  regular 
distortion  of  the  vertical  fringes  coincident  with  the  diagonal 
shocks  in  the  laser  shadowgram.  Density  profiles  across  the  tube 
at  station  7.5  diameters  behind  the  shock  confirm  that  the  shocks 
are  weak  (Figure  34). 

Viscous  effects  appear  to  be  confined  to  a  region  near 
the  wall  in  contradiction  to  the  main  assumption  in  one-dimensional 
computer  models.  The  shadowgrams  reveal  a  turbulent  region  only 
between  the  ribs.  The  interferogram  at  X/D  =  7.5  does  not  indicate 
a  large  cold  wall  influence  (Tm/Tw  =  11.8  for  M  =  7.5)  on  the 
flow  above  the  ribs.  Even  if  the  small  density  increase  near  the 
wall  in  Figure  34  is  interpreted  as  a  viscous  effect,  the  boundary 
layer  height  is  on  the  order  of  the  smooth  wall  boundary  layer 
thickness  (2.9  rib  heights)  predicted  by  the  method  of  Reference  25. 

The  inference  from  the  flow  visualization  results  is 
that  the  shock  flowfield  is  two-dimensional  in  the  vicinity  of 
the  main  shock,  and  that  two-dimensional  calculations  are  required 
to  understand  the  shock  attenuation  process.  The  influence  of 
viscosity  may  contribute  increased  attenuation  over  that  predicted 
by  inviscid  codes. 
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Figure  34,  Density  profiles  in  the  quasi-steady 
flow  (Pi  =  9,03  x  10“^  gm/cm^). 
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Wall  drag  measurements  were  performed  for  runs  60  and  61 
with  the  custom  built  load  cell  located  twenty  diameters  down¬ 
stream  of  the  entrance  to  the  MX  rib  test  section.  The  local  shock 
strengths  were  P2^  =  190  and  350  respectively.  Load  cell  signals 
(Figure  35)  are  characterized  by  a  sharp  initial  deflection  followed 
by  a  gradual  rise  to  peak  load  at  2  msec  after  shock  arrival.  The 
similarity  of  the  load  cell  and  pressure  signals  led  to  the  sus¬ 
picion  that  the  load  cell  was  sensitive  to  internal  pressure,  and 
this  was  confirmed  by  post  test  calibrations.  Vendor  estimates  of 
the  internal  pressure  sensitivity  indicated  that  at  the  shock 
arrival  the  internal  pressure  may  account  for  as  much  as  thirty 
percent  of  the  signal;  at  later  times  the  internal  pressure  could 
account  for  even  larger  fractions  of  the  signal.  Therefore,  only 
the  initial  signal  deflection  was  analyzed. 

Measured  drag  forces  were  divided  by  the  liner  wall  area 
and  a  dynamic  pressure  calculated  from 

q  =  2  P2  M2  (18) 

to  produce  an  effective  drag  coefficient  near  the  shock  wave 
(Table  3) .  The  steady  flow  drag  coefficient  according  to  equations 
1  and  2  is  0.032  which  is  a  factor  of  three  to  five  lower  than  the 
measured  values.  In  light  of  the  flow  visualization  observations, 
a  simple  wave  drag  model  was  formulated  which  at  least  produces 
higher  estimates.  Loads  on  the  MX  rib  liner  section  were  cal¬ 
culated  from  pressure  distribution  estimates  for  forward  and  rear¬ 
ward  facing  steps  neglecting  the  tangential  shearing  forces  on 
the  axial  portion  of  the  liner.  The  average  pressure  on  a  for¬ 
ward  facing  step  (Reference  21).  was  taken  as 

Pf  =  0.8  (M2  +  2)  P2  (19) 
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Figure  35.  MX  liner  load  histories 

(10  Klbs/div  X  1.0  msec/div), 
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and  the  average  pressure  on  a  rearward  facing  step  (Reference  22) 
was  taken  as 


P 


R 


(20) 


Both  correlations  are  for  data  for  stream  Mach  numbers  between  two 
and  four.  The  net  axial  force  on  the  rib  is  the  product  of  the 
pressure  difference  and  the  rib  frontal  area.  The  rib  force  is 
averaged  over  the  wall  area  per  rib  which  is  normalized  by  the 
dynamic  pressure  to  form  a  drag  coefficient 
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where  k/S  is  the  rib  height  to  spacing  ratio. 

The  few  experiments  conducted  here  are  not  conclusive. 
More  detailed  experiments  and  comparisons  with  the  models  are 
required  to  establish  an  effective  wall  drag. 


4 . 5  End  Wall  Reflected  Pressures 

The  end  wall  test  (Run  63)  was  run  at  an  initial  test 
section  pressure  of  0.1  atmosphere  and  an  initial  shock  strength 
at  the  rib  section  entrance  of  P2-l  =  80  to  minimize  the  loads  on 
the  shock  tube.  A  measurement  station  was  not  available  just  up¬ 
stream  of  the  end  wall,  so  the  shock  strength  incident  on  the 
wall  was  inferred  from  Run  64.  Initial  test  conditions  for  Run  63 
were  the  same  at,  for  Run  64,  and  pressure  measurements  at  four 
stations  upstream  of  the  end  wall  matched  well  for  the  two  runs 
(Figure  36) .  The  interpolated  value  of  the  incident  shock  strength 
is  P21  =  52,  and  the  theoretical  reflection  factor  for  a  plane 
shock  wave  of  this  strength  is  about  7.5  (References  15  and  17). 

End  wall  reflected  shock  waveforms  exhibited  an  irregular 
pattern  as  seen  in  the  plots  of  Figure  37  ^oscilloscope  traces  in 
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Figure  37,  End  wall  pressure  measurements. 


Appendix  B) .  The  three  pressure  transducers  measured  an  initial 
reflected  pressure  ratio  from  8.5  to  11.  The  duration  of  this 
unexpectedly  high  reflected  pressure  varied  from  20  ysec  at  the  outer 
transducer  (r  =  3R/4)  to  50  ysec  at  the  centerline  after  which  the 
pressures  decayed  to  near  the  expected  level. 

Pressure  fluctuations  apparently  not  associated  with  trans¬ 
ducer  ringing  (the  tranducers  have  a  natural  frequency  of  500  kHz) 
persist  throughout  the  records.  Interaction  of  the  reflected  shock 
and  the  nonuniform  region  behind  the  incident  shock  may  contribute 
to  the  pressure  fluctuations  for  about  200  ysec  after  shock  arrival. 

A  more  satisfactory  explanation  for  the  irregularities  and  seem¬ 
ingly  high  reflection  factors  is  apparent  from  the  optical  measure¬ 
ments  which  revealed  the  strongly  two-dimensional  nature  of  the  flow- 
field.  Secondary  shock  waves  created  by  the  ribs  cause  the  shock 
to  be  nonplanar  (small  differences  in  shock  arrival  time  were 
noted  on  the  end  wall)  and  create  pressure  variations  in  both 
the  radial  and  longitudinal  direction.  Thus,  pressures  mea¬ 
sured  on  the  rib  tops  may  not  be  representative  of  local  pressures 
across  the  tube. 


4 . 6  Shock  Attenuation  Comparison  with  Predictions 

Three  approaches  have  been  suggested  for  predicting 
shock  attenuation  in  a  ribbed  tube:  empirical  such  as  the  method 
due  to  Porzel,  one-dimensional  with  fully  developed  pipe  flow 
assumptions  and  two-dimensional  inviscid  calculations.  A  brief 
comparison  of  the  first  two  techniques  with  the  test  data  is 
presented  here.  Two-dimensional  calculations  are  being  performed 
by  Air  Force  Weapons  Laboratory  personnel. 

The  semi-empirical  prediction  method  of  Porzel  has  the 
simple  form 


(22) 


75 


for  strong  shocks.  Using  the  geometrical  rib  height  (k/D  = 

0.038)  Equation  (22)  does  not  predict  the  attenuation  observed 
in  the  present  tests  (e.g.,  Figure  38).  Although  an  expression  of 
the  form 


would  reproduce  the  test  results,  there  is  no  guarantee  that  the 
coefficient  is  constant  for  other  shock  strengths  or  other  rib 
designs.  Furthermore,  an  empirical  approach  such  as  this  ignores 
the  detailed  physics  of  the  attenuation  mechanism  and  does  not 
treat  the  entrance  region  effects  observed  in  the  preliminary 
tests . 

RIST  code  calculations  were  performed  for  the  tests  run  at 
initial  shock  strengths  of  P01  =  80  and  300  to  determine  if  one- 

^±i 

dimensional  shock  attenuation  models  based  on  steady  flow  assump¬ 
tions  were  adequate.  The  calculations  predicted  the  smooth  wall 
shock  attenuation  reasonably  well  (Figure  39)  ,  considering  that 
no  special  pains  were  taken  to  model  driver  conditions.  RIST 
also  did  an  acceptable  job  of  predicting  shock  attenuation  for  the 

lower  shock  strength  (P01  =  80) .  However,  RIST  code  calculations 

*-Li 

of  the  attenuation  for  P21  =  300  were  inadequate.  The  rapid  decay 
in  shock  strength  at  the  entrance  region  measured  in  the  initial 
tests  was  not  reproduced,  and  overall  shock  strengths  decayed 
about  twice  as  fast  as  predicted. 

No  attempt  at  improving  the  one-dimensional  code  was 
undertaken  in  light  of  the  entrance  effects  and  the  distinctly 
two-dimensional  flowfield  observed  in  the  optical  results.  Two 
dimensional  calculations  performed  by  Allen  Kuhl  for  the  TRW 
experiment  (Reference  20)  and  AFWL  personnel  for  the  present  ex¬ 
periment  need  to  be  studied  in  detail  before  further  modeling  is 
pursued. 
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Figure  38,  Comparison  of  Porzel  model  with  typical  data  set, 
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STREN6TH 


5.0 


SUMMARY 


A  shock  tube  test  program  was  conducted  at  the  NASA  Ames 
East  facility  to  provide  experimental  data  useful  for  MX  trench 
prototype  assessment.  The  experiments  considered  the  propagation 
of  strong  shock  waves  in  a  rigid  tube  with  transverse  internal  ribs 
In  addition  to  standard  time  of  arrival  and  pressure  measurements, 
custom  built  diagnostics  were  assembled  to  measure  wall  drag  and 
to  provide  an  optical  description  of  the  rib/f lowf ield  interactions 
Limited  comparisons  of  test  results  with  prediction  techniques 
were  made. 

5.1  Conclusions 


The  conclusions  reached  during  the  program  are: 

•  The  shock  tube  test  data  should  scale  to  the 
prototype  trench.  Experiment  scales  were  intro¬ 
duced  through  geometry  and  flow  Reynolds  number. 
The  experiment  was  performed  with  a  1:  39 
geometric  scaling  of  the  prototype  design.  The 
test  Reynolds  number  was  two  to  four  orders  of 
magnitude  lower  than  the  prototype.  How¬ 
ever,  no  discernable  effect  ot  flow  Reynolds 
number  on  shock  attenuation  or  on  pressure  wave¬ 
forms  was  seen  when  the  Reynolds  number  was 
varied  by  an  order  of  magnitude. 

•  Shock  wave  flowfields  are  two-dimensional 

in  the  vicinity  of  the  shock  front.  Shadow- 
grams  in  the  shock-front  vicinity 
revealed  a  complicated  but  symmetric  wave 
structure  composed  of  a  perturbed  main  shock 
and  reflected  shocks  originating  at  each  rib. 
Farther  behind  the  shock  (X/D  %  7)  the  flow- 
field  appeared  to  take  on  a  quasi-steady  char¬ 
acter  with  bow  shocks  emanating  from  each 
rib.  Interf erograms  at  X/D  =  7  produced  no 
evidence  of  fully  developed  pipe  flow. 

•  Shock  strength  attenuation  is  dependent  on 
initial  shock  strength  and  apparently  not 
strongly  dependent  on  rib  design.  High 
strength  shocks  (P21  >  300)  decayed  to  half 
their  smooth  wall  values  in  about  twenty 
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diameters;  lower  shock  strengths  decayed 
less  rapidly.  Entrance  effects  appeared 
to  play  a  strong  role  in  the  attenuation. 

An  alternate  rib  geometry  optimized  (with¬ 
in  design  constraints)  for  steady  flow  wall 
drag  produced  a  small  increase  in  shock 
attenuation. 

•  Two-dimensional  flcwfield  calculations 
are  required  to  understand  the  ribbed 
wall  shock  attenuation  process.  The 
empirical  formulation  of  Porzel  and  a 
one-dimensional  computer  code  model  do 
not  predict  the  observed  attenuations  for 
high  strength  shocks.  The  rapid  decrease 
in  shock  strengths  in  the  entrance  region 
appears  to  be  two-dimensional  and  must  be 
replicated  before  confident  prediction  of 
prototype  airblast  environments  can  be 
performed . 

•  Pressure  waveforms  in  the  ribbed  tube  are 
characterized  by  a  rise  in  pressure  behind 
the  shock  wave.  While  expansion  and  vent¬ 
ing  will  probably  minimize  this  effect  in 
the  MX  trench  prototype,  the  pressure  rise 
should  be  predicted  by  computer  codes  attempt¬ 
ing  to  calculate  the  experiments. 

5.2  Recommendations 


As  a  result  of  the  test  program  the  following  recommenda¬ 
tions  are  made  to  improve  the  understanding  of  shock  propagation 
in  tubes  with  transverse  ribs: 

•  Two-dimensional  flowfield  calculations  should 
be  performed  for  the  present  experiments  and 
closely  examined.  Special  emphasis  should  be 
placed  on  understanding  the  entrance  region 
phenomenon  to  determine  how  the  subsequent 
attenuation  is  influenced  and  whether  it  is 
test  unique.  The  calculations  should  be 
compared  with  measured  shock  time  of  arrival, 
shock  strength,  pressure  waveforms,  flow 
structure  and  density  distribution. 

•  The  possibility  of  contructing  simplified 
models  of  the  shock  attenuation  process 
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should  be  examined.  One  dimensional  com¬ 
puter  codes  provide  a  cost  effective  design 
tool  for  computing  shock  propagation  if 
adequate  models  can  be  developed. 

Design  tools  should  be  verified  with  se¬ 
lected  experiments.  When  a  shock  propa¬ 
gation  methodology  has  been  developed,  the 
model  can  be  exercised  over  a  range  of 
shock  tube  parameters  to  define  ex¬ 
periment  conditions  which  would  test  the 
universality  of  the  method. 

A  shock  attenuation  wall  treatment  tor 
hardened  base  entrainment  systems  should  be 
developed  and  tested  based  on  the  knowledge 
gained.  An  effective  shock  attenuation 
system  could  relieve  the  design  requirement 
on  blast  valves  or  doors. 
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APPENDIX  A 


Time  of  Arrival  Data 

Shock  wave  time  of  arrival  (TOA)  data  from  the  MSTIMS 
system  (smooth  tube)  and  from  oscilloscope  traces  of  pressure 
signals  (ribbed  tubes)  are  catalogued  in  Tables  A-l  to  A-4. 

Time  is  measured  from  180  ysec  after  peak  current  of  the  arc  dis¬ 
charge  and  is  accurate  to  5  usee.  Station  location  is  measured 
from  the  diaphragm.  The  order  of  the  presentation  is  consistent 
with  Table  1  in  the  text. 


A-l 


o 

VO 


Table  A-4.  Time  of  arrival  data  for  Series  A, 
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Roughness 


APPENDIX  B 


Pressure  Data 

Signals  from  the  pressure  transducers  were  displayed  on 
Tektronix  oscilloscopes  and  recorded  on  Polaroid  film.  Polaroid 
reproductions  for  the  final  runs  are  presented  in  Figures  B-l 
through  B-15;  corresponding  identification  and  sensitivity  in¬ 
formation  is  catalogued  in  Tables  B-l  through  B-15.  Shock  wave 
pressure  ratio  measurements  determined  from  the  initial  signal 
deflection  peak  (the  midpoint  of  the  trace  was  used)  are  pre¬ 
sented  in  Tables  B-16  through  B-19  consistent  with  Table  1  of 
the  text.  The  uncertainty  in  reading  the  signal  from  the  photo¬ 
graphs  is  taken  to  be  the  ratio  of  half  the  beam  width  to  the 
initial  beam  deflection  which  is  generally  less  than  10  percent. 
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B-2.  OSCILLISCOPE  SETTINGS  FOR  RUN  7 
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Table  B-5.  Oscilliscope  settings  for  Run  52. 


Table  B-6.  Oscilliscope  settings  for  Run  54. 
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Table  B-7.  Oscilliscope  settings  for  Run  56. 


oooo  oooo 
oooo  oooo 
in  in  in  in  in  in  in  in 


m  in  | 


oooo  oooo  in  in  I  in  in  I  in  in 

rHrHrHf-i  iH<N<N<N  O 


id  h  id  r^vor^o  vovo  vo  r-  vo  r* 

oooo  o^rch^r  oo  oo  oo  ^  o\ 

rHiHrHr— i  H  lH  ^  ID  iHrH  r-H  t— I  H  H  in 


H  CN  H  CN  H  N  H  (N 

P  P  J  PPPJ 


rl  M  CO  ^  If)  VO  CM  rH  CO  CO 

(M  Q  K  «  «  «  K  «  Q«  Cti«  >« 


a)  a>  a)  <u  <u  <u  <d  <u  q>  a>  <u  <d  I  a>  <d  a>  <u  <u 


MM  MM 
3  3  3  3 


MMMM  MMMM  MM  MM  MM  MM 

3333  3333  33  33  33  33 

CO  01  CO  03  03  03  01  03  COCA  COO)  COCO  cool 

cocococo  co  co  coco  coco  coco  coco  coco 

cucucuai  <u  a)  o)  a)  ojcu  <u  <d  a><u  a)  <u 

MMMM  MMMM  MM  MM  MM  MM 

CM  0-.  CM  CM  ft  ft  Pi  ft  ft  CM  CM  CW  CM  P-t  pH 


MMMM 
3  3  3  3 

CO  CO  CO  CO 

CO  CO  CO  CO 

a)  a>  a)  a) 

MMMM 


coco  coca  coco  coco 

aja)  q)  at  a>  a)  aj  q) 


56- 


B-16 


Table  B-8.  Oscilliscope  settings  for  Run  57, 
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Table  B-9.  Oscilliscope  settings  for  Run  58. 
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igure  B-lu.  Pressure  signals  for  Kun 


Table  B-10.  Oscilli scope  settings  for  Run  59. 
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Table  B-13.  Oscilli scope  settings  for  Run  60, 
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Tarle  B-16.  Shock  front  pressure  ratio  data  for  Series  1. 


APPENDIX  C 


Optical  Measurements 

This  Appendix  presents  a  discussion  of  an  optical  in¬ 
vestigation  into  the  complex  flow  about  the  ribs  in  a  model  of 
the  MX  trench.  Using  holographic  techniques,  a  pictoral  de¬ 
scription  of  the  wave  and  flow  structure  has  been  compiled  for 
various  times  after  shock  arrival.  Photographs  of  the  holo¬ 
graphic  images  obtained  using  a  variety  of  methods  are  included 
along  with  a  discussion  of  interferometric  density  measurements. 

An  SAI  custom  built  laser  holography  system  was  used  to 
record  several  types  of  flowfield  images  including  laser  shadow- 
graphy,  conventional  inline  and  off-axis  holography,  and  holo¬ 
graphic  interferometry.  For  this  system,  the  source  of  coherent 
illumination  was  an  Nd:YAG  frequency  doubled  laser  emitting  30 
millijoules  of  5230  A  light  in  20  nanosecond  pulses. 

The  optical  system  is  shown  schematically  in  Figure  C-l. 
The  holographic  camera  has  the  capability  of  focusing  the  image 
of  the  center  of  the  test  section  into  the  hologram  plane  when 
laser  shadowgraphy  is  employed;  or  optionally  focusing  the 
image  into  a  plane  upstream  from  the  holograph  plane  when  off- 
axis  holography  is  employed.  An  optional  off-axis  reference 
beam  provided  the  capabilities  to  perform  off-axis  holography. 
The  laser  was  externally  spatially  filtered  to  provide  a  smooth 
beam,  and  the  entire  laser  train  was  rail  mounted  for  stability. 

For  laser  shadowgraphy,  a  pulse  laser  provided  an 
adequately  short  (20  nanosecond)  exposure.  Although  it  is  not 
essential  to  focus  an  image  of  the  central  area  of  the  flow 
into  the  film  plane,  auxiliary  optics  were  used  to  obtain  in¬ 
focus  images  (Figures  C-2  thru  C-4) .  These  three  shadowgrams 
show  the  shock  as  it  appears  within  the  window  at  three  dif¬ 
ferent  times  during  the  event.  The  corrugated  shock  tube  walls 
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Figure  C-l.  Holography  system. 
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Figure  C-2.  Laser  shadowgraphy  (Run  38,  t  =  12  ^sec). 


evident  in  the  upper  and  lower  portions  of  the  photograph  give 
rise  to  a  complicated,  symmetrical,  regular,  and  relatively  tur¬ 
bulence  free  shock  structure.  The  photographs  demonstrate  that 
certain  gross  features  of  the  shock  can  be  obtained  with  simple 
laser  shadowgraphy. 

Laser  shadowgraphy  can  be  considered  as  a  form  of  in-line 
holography.  In-line  holography  uses  a  clean  reference  beam  which 
is  caused  to  interfere  with  the  object  transmitted  light,  and 
whose  propagation  direction  is  identical  to  the  major  portion  of 
the  object  wave  (i.e.,  no  offset  angle).  In  the  shock  tube  con¬ 
figuration,  the  reference  wave  may  be  regarded  as  that  portion 
of  the  illumination  wave  undiffracted  or  otherwise  unaltered  by 
the  object.  If  this  wave  is  of  significantly  greater  amplitude 
at  the  film  plane  than  the  object  modified  portion  of  the 
illumination  (a  necessity  for  nonlinear  recording) ,  a  useful 
hologram  is  constructed.  An  object  field  containing  a  dense 
cloud  of  particles  will  not  do,  nor  will  one  in  which  the  tur¬ 
bulence  (through  which  the  illumination  beam  must  pass) is  too 
great.  Depending  upon  the  ultimate  measurement  to  be  made  from 
the  hologram,  the  results  (Figures  C-2  through  C-4)  are  either 
acceptable  or  unacceptable  in  line  holograms.  There  are  few 
particles  in  the  field,  but  there  is  some  turbulence  near  the 
tube  walls.  Upon  closer  examination,  one  would  find  that  the 
illumination  wave  is  diffracted  in  a  complex  manner  by  the 
shock  wave.  Another  disadvantage  in  using  in-line  holograms 
is  that  a  conjugate  image  overlaps  the  desired  image  when  the 
hologram  is  reconstructued,  thus  often  contributing  to  confusion 
in  data  interpretation. 

The  limitations  of  in-line  holography  prevent  one  from 
using  the  technique  to  measure  the  density  profiles  within  the 
shock.  Such  measurements  require  the  recording  of  the  phase  of 
the  light  for  each  point  across  the  film  plane,  so  that  inter¬ 
ferometry  may  be  implemented. 


Good  quality  holography  is  obtained  only  when  the 
reference  wave  is  clean.  Holograms  constructed  with  an  off- 
axis  holography  system  (Figure  C-l)  give  rise  to  the  recon¬ 
structed  images  shown  in  Figures  C-5  and  C-6.  These  photographs 
are  interesting  in  that  they  show  how  well  a  duplication  of 
the  previous  shock  phenomena  is  achieved.  They  also  demonstrate 
that  the  photographic  plates  used  to  record  the  events  for  Fig¬ 
ures  C-2  through  C-4,  were  somewhat  nonlinear  when  compared  with 
the  film  results  of  Figures  C-5  and  C-6  which  lack  the  contrasting 
dark  and  light  lines. 

Holograms  constructed  with  the  present  system  contain 
records  of  both  the  amplitude  and  phase  of  the  object  wave  (the 
essence  of  holography) ,  and  the  reconstructed  image  waves  can 
be  allowed  to  interfere  with  a  reference  wave  to  create  inter¬ 
ference  fringes.  The  fringe  structure  contains  information  from 
which  can  be  extracted  details  about  the  density  for  all  points 
within  the  shock  wave.  Rather  than  attempt  to  carry  out  this 
procedure,  a  second  method  of  interferometry  was  implemented 
which  has  both  advantages  and  disadvantages  over  this  post 
event  method.  This  second  method  requires  that  one  record  two 
exposures  on  the  single  photographic  plate.  One  of  the  exposures 
represents  the  reference,  and  thus,  the  interferometric  structure 
is  contained  within  the  single  holographic  image.  An  obvious 
disadvantage  of  the  method  is  that  the  interferometric  fringe 
structure  is  frozen  and  cannot  be  altered  during  the  reconstruc¬ 
tion  and  interrogation  steps.  An  advantage  is  that  the  effects 
of  the  distortions  and  aberrations  caused  by  the  optics  in  the 
construction  apparatus  are  automatically  cancelled  from  the 
resulting  fringe  structure. 

An  example  of  interferometry  with  zero  offset  in  which 
the  first  exposure  occurred  prior  to  shock  propagation  and  the 
second  exposure  occurred  during  shock  propagation  is  shown  in 


in  Figure  C-7.  The  image  is  of  an  event  which  is  similar  to 
those  of  Figures  C-2  and  C-5.  The  complicated  fringe  structure 
is  of  such  a  nature  that  interpretation  is  difficult,  if  not 
impossible.  It  is  not  clear  whether  adjacent  fringes  represent 
regions  of  falling  or  rising  density.  A  method  which  overcomes 
this  problem  is  the  following.  Between  the  first  and  second 
exposures,  the  illumination  wave  (or  alternatively  the  hologram 
construction  reference  wave)  direction  is  slightly  altered  so 
that  a  residual  vertical  fringe  structure  overlays  the  holo¬ 
graphic  image  and  reduces  the  ambiguities.  The  degree  of  direc¬ 
tional  change  determines  the  residual  fringe  spacing;  too  great 
a  tilt  causes  fringes  which  are  too  fine  to  be  useful  while  not 
enough  tilt  contributes  little  to  reducing  ambiguities.  This 
procedure  results  in  offset  fringe  interferograms  with  vertical 
fringes  (Figures  C-8  through  C-ll) .  Three  of  these  (Figures  C-8 
through  C-10)  were  recorded  at  approximately  the  same  relative 
time  as  the  laser  shadowgrams. 

Density  profiles  are  derived  from  the  fringe  structure 
by  relating  fringe  shift  to  index  of  refraction  changes (which 
for  gases  is  related  to  the  density  change).  From  the  Lorentz- 
Lorenz  formula  (Reference  C-l)  the  polarizability  a  of  the  gas 
within  the  shock  volume  is 


a 


(C-l) 


in  which  the  number  of  molecules  per  unit  volume  N  and  the  index 
of  refraction  n  are  introduced.  The  index  of  refraction  of  a 
gas  is  close  to  unity  even  at  high  pressures;  therefore.  Equa¬ 
tion  C-l  reduces  to 


n  -  1  =  cp 


(C-2) 
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Figure  C-9.  Douri  f  piii  sf  i  asfr  holography  offsft 
fringf  (Run  44,  t  =  73/iSEc). 
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gure  C-ll.  Double  pulse  laser  holography  offset 
fringe  (Run  46 ,  t  =  341^sec). 


in  which  the  density  p  and  a  proportionality  constant  c  with 
units  of  inverse  density  are  introduced. 

A  nondimensional  constant  can  be  introduced  by  defining 
a  standard  density.  This  number  is  the  familiar  Gladstone  Dale 
constant  $.  Equation  C-2  is  then 

n  -  1  -  e  g-  (C- 

Ms 

-4  -  3  3 

where  3  =  2.95  x  10  and  p  =  1.29  x  10  gm/cm  . 

It  remains  however  to  relate  the  index  of  refraction  to 
the  interferometric  fringe  spacing.  In  the  absence  of  an  offset 
or  tilting  between  the  two  exposures  (or  alternatively  when  cal¬ 
culating  variations  in  a  direction  along  a  linear  fringe  direc¬ 
tion,  i.e.,  top  to  bottom  in  Figures  C-8  through  C-ll) ,  the  index 
of  refraction  variation  between  adjacent  fringes  An  is 

An  =  p  (C- 

where  the  test  section  width  D  (the  mean  optical  path  through 
the  test  section)  and  the  wavelength  of  the  laser  light  X  are 
introduced. 

The  density  change  between  adjacent  fringes  is  related 
to  the  index  of  refraction  change  by 

Ap  =  3H  *  An  “  (m)  ps  (c" 

For  these  tests,  the  actual  interferometer  sensitivity  is 

1  ps  -53 

Ap  *  21 A  fringe  =  4*71  x  10  gm/cm  -fringe  (C- 

Density  profiles  were  obtained  by  counting  fringes  outwards  from 
the  centerline  of  the  test  section.  The  absolute  centerline 


density  was  estimated  from  the  shock  speed  and  the  initial  den- 

-3  3 

sity  to  be  9.03  x  10  gm/cm  ,  and  the  density  profiles  are 
plotted  relative  to  the  centerline  value  (Section  4.5). 
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